Pure Fulde-Ferrell-Larkin-Ovchinnikov state in optical lattices 
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We study the phase diagram of one-dimensional, two-component attractive fermions on optical lattices with 
an olf-diagonal confinement. We identify in this system a pure Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
with spatially modulated superfluidity and polarization in a large parameter window, which provides an easier 
playground to detect the FFLO state experimentally. The clear signature of FFLO state is analyzed in the 
ground state properties, the pair correlation function, and the magnetic structure factor, using the density-matrix 
renormalization group method. 
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Introduction. — Experimental developments in cold atomic 
gases provide a promising testbed to realize the exotic quan- 
tum phases predicted in condensed matter systems (for details 
see Ref. [1] for a recent review). Among these phases, there 
is an important state first predicted by Fulde, Ferrell, Larkin, 
and Ovchinnikov (FFLO) in 1964 for superconductors in high 
magnetic fields, characterized by cooper pairs with nonzero 
center-of-mass momentum and a spatially inhomogeneous or- 
der parameter.^ 

Though the FFLO state has not been explicitly confirmed 
experimentally, its realization in ultracold atomic systems is 
becoming promising due to the great adjustability for almost 
all of the system parameters, for example, the spin imbal- 
ance realized by the magnetic fields in solid state physics can 
be easily simulated using the mismatched different particle 
species by radio-frequency fields.^ While the FFLO state is 
difficult to detect in three-dimensional ultracold atomic sys- 
tems due to the very small parameter window in free spaced 
the possibility of direct realization in the one-dimensional 
(ID) case increases due to the major part of the phase dia- 
gram being of the character of the FFLO correlation, based on 
analytical,^ '' numerical,^ and experimental*^ research. 

Experimentally*^ the density profiles of a spin imbalanced 
mixture of ultracold ^Li atoms trapped in an array of ID tubes 
are measured with the standard absorption imaging. The den- 
sity profiles detected meet the criteria of the FFLO state pre- 
dicted analytically.^ In a ID system of harmonic confinement, 
a partially polarized core is surrounded by wings of either 
a fully paired or a fully polarized Fermi gas depending on 
the controlled imbalanced population in its two spin states. 
The partially polarized core in the density profile is proved of 
FFLO characteristics by calculating the pair momentum.* But 
in order to demonstrate directly that this pair state is of FFLO 
type, one needs to precisely measure the physical quantities 
peculiar to FFLO states, for example, the condensed fraction 
using the time-of-flight (TOF) measurement and the quanti- 
ties other than density profile revealing the finite momentum 
of the pairing order parameter, like the correlations in the shot 
noise of fermion absorption images in TOF^ '^"^ In a ID op- 
tical lattice under the presence of a harmonic confining po- 
tential, Korolyuk et al. proposed that the change of the dou- 
ble occupancy due to the lattice depth modulation provides 
a clear evidence." The collective modes produced by the re- 
sponse of the ground state to time-dependent potentials give 



another possible signature of the FFLO phase. '^ '-^ Very re- 
cently, Kajala et al. proved that the TOF expansion manifests 
a direct signature of the FFLO state in the ID system in which 
the expansion velocity of the unpaired particles matches the 
expected FFLO momentum. 

In cold atom experiments, the necessary confining poten- 
tial, which is usually a single-particle potential called diago- 
nal confinement, leads to the emergence of coexisting states'— 
and inevitably makes the detected signals complicated. In 
the repulsive case the diagonal confinement leads to the ab- 
sence of true incompressible Mott insulating phases and a real 
superfluid-to-Mott transition."" Though Scarola et al. pointed 
out that edge effects in a trapped Fermi-Hubbard system can 
be filtered out with a direct probe independent of inhomo- 
geneity'^ by relating the core compressibility to changes in 
the double occupancy, people hope to realize directly a pure 
Mott insulating or other exotic phases in the trapped system 
by making use of state-of-the-art cold atomic techniques. In 
order to achieve such a goal, an off-diagonal confinement 
method was pointed out by Rousseau et al.r^il^ where site- 
dependent hopping integrals are designed with amplitude de- 
creasing while moving away from the center of the lattice and 
vanishing towards the boundaries. 

In this paper, we show that for the system in the off- 
diagonal confinement of population imbalance, a pure FFLO 
state is realized at a large range of system parameters, and 
thus, in this case, the common multi-phase structure in the di- 
agonal confinement system is avoided, which facilitates the 
detection of such an elusive state in the cold atomic system. 
The clear signature of FFLO state is analyzed through calcu- 
lating the pair correlation function and the magnetic structure 
factor, using density-matrix renormalization group (DMRG) 
method. 

The model and phase diagram. — We consider fermions 
confined to a ID optical lattice with L sites and lattice con- 
stant a — 1. The Hamiltonian is given by 
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where cl^J is a fermionic annihilation (creation) operator act- 
ing on site i with pseudospin a representing the hyperfine 
states and the sum over (i, j) is restricted to the nearest neigh- 
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hours. The spatially varying tunneling is of inverted parahola 
designed as 

Uj = t{i+j + 1)(2L -I- J- (2) 

which vanish at the edges avoiding the escape of the particles 
from the edges, namely, i-i,o = th-i.L ~ 0. The properties 
of the FFLO state in a ID optical lattice of mass imbalance 
have been studied intensively,'^^' where in Eq. (1) is re- 
placed by tcr- They found that the typical FFLO modulation 
is preserved in the gapped phase but not the quasi-long-range 
order'^ and if the system is partially polarized, the ground- 
state FFLO state is more stable in the presence of mass imbal- 
^1 

ance.- 

We study the case of attractive on-site interactions (U < 0) 
with population imbalance. In the following, we take u — U/t 
with t as the units of energy. The number of fermions of 
each species is = '^i{cfc,Ci,a) with N ~ N-^ + 
and the spin a =t, I. The filling factor, the magnetization, 
and the polarization are n — N/L, Sz = {N^ — Ni)/2, and 
p — {N-i-~Ni)/N, respectively. We work in the canonical en- 
semble with fixed total number of particles. In the following 
we choose a typical experimentally accessible optical lattice 
of L = 70 sites and assume iV^ > without loss of gener- 
ality. 

Experimentally, this off-diagonal confinement system can 
be realized by a holographic method with a mask under the 
present-day ultracold atomic techniquesji^ Technically, the 
off-diagonal confinement system is somehow equivalent to the 
smooth boundary conditions introduced by Vekic and White 
aiming at reducing finite-size effects and extrapolating to the 
thermodynamic limit on relatively small systems, and is 
similar to the sine-squared or sinusoidal deformation applied 
to ID quantum Hamiltonians.^-' In our study we discuss a par- 
tially polarized system of inverted parabola tunneling.^"' 

The ground-state phase diagram is shown in Fig. [T] The 
typical density distributions in different ranges are displayed 
as insets. We find that a partially polarized phase (0 < 
ni^-, nil < 1) exists in the core of the system at finite polar- 
ization. Among all the different phases, we identify an impor- 
tant pure FFLO superfluid phase which is partially polarized 
in the whole space of the system in part I, for Sz £ [8, 18] 
and n G [0.93, 1.1]. We would like to point out that in the 
diagonal confinement system the pure FFLO state happens 
only at a certain critical polarization. Other phases identified 
are, a fully polarized phase (0 < fii^ < l,nii = 0, FP), 
an unpolarized BCS phase (n^^ = ^ 0, BCS), a pure 
Fock state of localized pairs, namely a band-insulating state, 
at high densities (n^^ = n^^ = 1, FS), and fully polarized 
ferromagnetic insulating phase for one species (riif — 1) and 
partially filled {riii ^ 0, Fla) or empty (rij^ = 0, Fib) for 
another. All the regions in Fig.[T]except part I are composed 
of coexisting phases, either of two-phase structure: FFLO-FP 
(Ila, lie), FFLO-FIa (Ild), and FFLO-BCS (Ilia, IIIc), FFLO- 
Flb (IV), or of three-phase structure: FFLO-FI-FS (lib) and 
FFLO-BCS-FS (Illb). The relation between the parts of Ila 
and lib, lie and lid, and Ilia and Illb, can be understood by 
the particle-hole symmetry, da (— l)'c^, mapping the 
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FIG. 1: (Color online) The ground-state phase diagram of the off- 
diagonal confinement system in the Sz-n plane for it = — 4. The 
pure FFLO superfluid phase (I) locates in the parameter range of 
8 < 5*2 < 18 and 0.93 < n < 1.1. Other composite phases we 
labeled are FFLO-FP (Ila, lie), FFLO-FI-FS (lib), FFLO-FIa (lid), 
FFLO-BCS (Ilia, IIIc), FFLO-BCS-FS (Illb), and FFLO-FIb (IV), 
respectively (details see text). The density distributions for spin-up 
and spin-down particles in the different regions are shown versus i 
for typical situations as insets. The arrows indicate the direction of 
the different phases while increasing the magnetization Sz for fixed 
fermions of A*' = 70 and 80. 



Hamiltonian (1) onto itself. As a result, if the particle number 
in part Ila and part lib satisfies + N\ = ^1 iV^ and 
N°- + = 2i, the densities in the corresponding parts are 
related by nf^ = 1 — n^^. Here tr — —a. 

Pairing correlations. — Let us now come to the discussion 
of the pairing properties. We first compute the pair correlation 
function, which is defined as, 

G(i,j-) = (A|A,) (3) 

with Ai — Ci^-Cii- Its Fourier transform is given by 

Gpair(fc) = ^EG(z,jV'^(-^). (4) 

hj 

Feiguin and Heidrich-Meisner confirmed that for the attrac- 
tive Hubbard model in a harmonic confinement the momen- 
tum pair distribution has peaks at the mismatch of the Fermi 
surfaces. This quantity is displayed in Fig. |2] for varying 
magnetization from 5^ = to 30. 

For the unpolarized case, the ground state is a BCS state 
characterized by a sharp peak centered at momentum fc = 
in the pair momentum distribution Gpair(fc)- For the pop- 
ulation imbalanced Fermi system of two spin species, two 
different Fermi surfaces exist and the ground state is a ID 
FFLO state with a finite center-of-mass momentum (fc ^ 0), 
which can be understood as an order parameter of the sys- 
tem. In the homogeneous system, the momentum of the FFLO 



- S^=1 1 
-S =15, 



-20 20 



_S^=2 
- S =6 



-20 20 



^.•11 

J'" Pure FFLO (I) 




-20 20 




eff 



FIG. 3: (Color online) Upper panel: the pair momentum distribution 
for different interaction strength. The parameters used here are Sz = 
15 and A'^ = 70. The one in the lower-left magnifies the region 
around the position of the peaks. Lower panel: the corresponding 
density difference An^ = riif — in the real space. The notations 
are the same as those in upper panel. 



FIG. 2: (Color online) Upper panel: the pair momentum distribu- 
tions for the magnetization varying from Sz = to 30 (from top to 
bottom) in a system of = 80 fermions. The inset shows that the 
momentum of the FFLO state fcppLO, at which the Gpair is peaked, 
as a function of the effective magnetization at the core S^^/Leg . 
With the increase of Sz, the phases of the system change from 
Illb ^ lib ^ Ild ^ IV for iV = 80 and from IIIc ^ I ^ IV 
for N = 70 (in the lower panel), which is indicated in Fig.[T]by ar- 
rows. The dashed line represents fcppLO = Afcp. Lower panel: the 
position of the FFLO peak as a function of /icff in a system of 
N — 70 fermions. fcppLO in the pure FFLO region is denoted by the 
full circle. In the inset, the profiles of the populations difference An^ 
are shown for two different magnetizations Sz in the corresponding 
phase region. The crosses with Sz values are pertinent to the insets. 
See the text for a discussion of the relations between the linearity of 
fcppLO and Arii in different phases. The interaction strength used 
here isu — —4. 



state fcpFLO^ at which the C?pair is peaked, is, /cfflo = 
/cpt ~ = 7r(n-f- — ni) = 2ttSz/L — nnp ^ with 
kpa = T^ria = irNa/L. In the inhomogeneous system, we 
define Afcp = / Lcs, where is obtained by inte- 

grating {SD over the effective partially polarized FFLO re- 
gion Leff with S], defined as 5* — (riif — nii)/2}^ We notice 
that, in Fig. |2] the value of fepFLO increases with Sz as pre- 
dicted by the FFLO picture. But in the inhomogeneous system 



under the off-diagonal confinement, where the hopping ampli- 
tude tij is spatially varying, the relation between /cfflo ™d 
Afcp is not as simple as that under the diagonal confinement. 
For N = 80, in the cases where the density differences in the 
bulk become inhomogeneous, fcppLo deviates slightly from 
Afcp. For the system of iV = 70 (lower panel), the phases of 
the system undergoes IIIc ^ I IV while increasing Sz- 
Starting from Sz = 13, fcppLO greatly deviates from Afcp at 
which Arii = — n^^ in the bulk becomes largely inho- 
mogeneous. The FFLO peak around fcppLo comes from the 
FFLO state in the bulk of the density, while the central peaks 
(of oscillating character) come from the paired BCS wings. 

For the pure FFLO states, we show in Fig. [3] the effects of 
the interactions on the pair momentum distributions. We find 
the optimal interaction strength (u « —15) where the peak of 
Gpair reaches maximum. The position of the peak matches 
exactly the linearity relation of fcppLO = 2ttSz/ L w 1.346 at 
u « —10, where An, is almost homogeneous, in the whole 
range of the system except at the edges. 

Magnetic structure factor. — We further reveal fingerprints 
of the FFLO state by analyzing the magnetic structure factor 
(MSF), which is the Fourier transform of spin-spin correla- 
tions, 

- ^ Y.^{SlSi) ~ {Sl){Si))e^'^''-^^ . (5) 
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FIG. 4: (Color online) Upper panel: magnetic structure factor as a 
function of k/n with fixed number of fermions N=70 and interaction 
strength it = — 10 for different magnetization Sz (from top Sz ~ 
to bottom Sz = 30). Middle panel: the corresponding momentum 
derivative of the magnetic structure factor. Lower panel: correspond- 
ing density imbalance distributions, 2(ni^ — riii). The dashed lines 
mark the matching positions of the kink in the structure factor, the 
corresponding jump in the derivative of the magnetic structure fac- 
tor, and the central density imbalance. 

Spin-spin correlations can be detected in a non-destructive 
way via spatially resolved quantum polarization spectroscopy. 
Roscilde et al. showed analytically and numerically that 
2fcFFLO-modulation of spin-spin correlations in the MSF is a 
direct fingerprint of pairing in an imbalance system and serves 
as a good quantity in detecting the FFLO state. 

We display the MSF in Fig. |4] for TV = 70 and w = -10 



with different polarization values. Comparing to the unpo- 
larized case, there is an obvious kink in the MSF with finite 
Sz ■ To have a clear signature of this kink, we show the mo- 
mentum derivative of the MSF which exhibits a marked jump, 
where the density imbalance distributions in the bulk, An^, 
reveal a constant structure. The match of the location of the 
kink with the bulk region of the density deteriorates when the 
density difference becomes inhomogeneous, like the cases in 
the lower panel of Fig. |2] when Sz > 13. We notice that the 
density difference in the lower panel is obviously oscillatory 
for small polarizations. The wavelength of the oscillations is 
given by A = 27r/|fcF^ — fepj.], which gives a confirmation 
that the FFLO state is inhomogeneous and exhibits real-space 
modulationsi2i 

Conclusions. — We have considered a ID optical lattice 
of off-diagonal confinement, modeled by a single-band Hub- 
bard model with spatially varying hopping amplitudes. We 
showed that with population imbalance the dominant pairing 
mechanism in the ground state is FFLO characterized by a fi- 
nite center-of-mass momentum. We predicted based on the 
Sz-n phase diagram that a pure FFLO state exists at a large 
parameter range, greatly different from the diagonal confine- 
ment system where only a critical polarization can achieve 
that. We further revealed that the partially polarized den- 
sities are the results of the finite center-of-mass momentum 
pairing through calculating the pair correlation function and 
the spin-spin correlation. The flexibility in ultracold atomic 
experiment could be used in designing the spatially varying 
hopping amplitude, and thus, provides an ideal playground 
for realizing such a FFLO superfluid state. Furthermore, the 
existence of the pure FFLO phase could facilitate an unam- 
biguous experimental observation of FFLO pairing in atomic 
systems free of the influence from other quantum phases at 
edges. The MSF which can be detected by the quantum polar- 
ization spectroscopy signal also provides a direct signature of 
finite center-of-mass momentum FFLO pairing. 
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